Sheet-current distribution in a dc SQUID washer probed by vortices 
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We present a novel method, based on vortex imaging by low-temperature scanning electron mi- 
croscopy (LTSEM), to directly image the sheet-current distribution in YBa2Cu3C>7 dc SQUID wash- 
ers. We show that the LTSEM vortex signals are simply related to the scalar stream function de- 
scribing the vortex-free circulating sheet-current distribution J. Unlike previous inversion methods 
that infer the current distribution from the measured magnetic field, our method uses pinned vor- 
tices as local detectors for J. Our experimental results are in very good agreement with numerical 
calculations of J. 
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Spatially resolved techniques can provide important 
insight into current flow, arrangement of vortices, flux 
pinning, and noise in superconductors and their mutual 
interactions. So far there has been only one method of 
imaging the current distribution in superconductors: The 
magnetic field distribution on top of a superconducting 
thin film is measured, e.g. by magneto-optics, from which 
the current distribution can then be calculated by invert- 
ing the Biot-Savart law[l|. 

In this paper we present a novel method to directly 
image the sheet-current distribution in a YBa2Cu3C>7 
thin film. We use low-temperature scanning electron 
microscopy (LTSEMjIlJ! i, |3 to image vortices in 
dc SQUID washers ja 0- Most techniques for vor- 
tex imaging, such as Lorentz microscopy |8j, scanning 
SQUID microscopy 01, fiof. scanning Hall microscopv|ll| 
or magneto-optics [T3 rely on the detection of the stray 
magnetic field produced in close proximity to a vortex. 
In contrast, vortex imaging by LTSEM is different from 
those techniques, as it is based on the electron-beam- 
induced apparent displacement of a vortex, pinned at 
position r in the (x, y)-plane of a SQUID washer, which 
is detected as a change of stray magnetic flux $(r) cou- 
pled to the SQUID. Hence, the contrast of the LTSEM 
vortex signals directly senses V$(r). Recently, Clem and 
Brandt |13) have shown that $(r) is proportional to the 
scalar stream function G(r) that describes the circulat- 
ing sheet-current density J{r) flowing in the vortex-free 
case at position r in the SQUID washer. In this paper we 
show that this relationship allows us to use the vortices 
as local detectors for J(r): At each position a vortex has 
been imaged, we can directly determine J(r) without 
complicated calculations. 

In our experiments, we investigated several dc SQUID 
washers [see Fig. da)] fabricated from epitaxially grown 
<i=80nm thick c-axis oriented YBa2Cu307 (YBCO) thin 
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films. We will present an analysis of LTSEM data ob- 
tained from one representative device with washer size 
120 /im x 305 /im, with a 100 /im long and 4 /im wide 
slit. The 1 /xm wide Josephson junctions are formed by a 
24° symmetric grain boundary in the underlying SrTi03 
substrate. For imaging by LTSEM, the YBCO SQUIDs 
are mounted on a magnetically shielded, liquid nitrogen 
cooled cryostage of an SEM [14j and read out by a stan- 
dard flux-locked loop (FLL) with 3.125 kHz bias-current 
reversal to eliminate 1// noise due to fluctuations in the 
critical current I c of the Josephson junctions. All LT- 
SEM images are obtained at T = 77 K after cooling the 
SQUID through T c in a static magnetic field Bq. 

For the spatially resolved measurements, the electron 
beam is used as a local perturbation that induces an in- 
crease in temperature 5T(x — xq, y — yo) on the sample 
surface (in the x, y-plane) centered on the beam spot po- 
sition (xo, yo)- At T=77K the length scale for the spatial 
decay of the thermal perturbation is set by the beam- 
electron range R « 0.5 /im for a typical beam voltage 
Vf, = 10 kV 4j. This gives a maximum increase in beam- 
induced temperature AT of a few K at (xo, yo) for typical 
beam currents lb of a few nA. So-called 6$(xo,yo) im- 
ages are obtained by recording the e-beam-induced flux 
change (5<f> in the SQUID as a function of the e-beam co- 
ordinates (xo,yo)- To improve the signal to noise ratio, 
we use a beam-blanking unit operating at typically 5 kHz, 
and the FLL output signal, i.e., the e-beam-induced flux 
change in the SQUID, is lock-in detected. 

The mechanism of imaging of pinned vortices can be 
briefly described as follows: The e-beam-induced local 
increase in temperature produces a local increase in the 
Pearl length [15| A = \\/d, where Xl is the temperature- 
dependent London penetration depth. Hence, the screen- 
ing currents circulating around a pinned vortex are spa- 
tially extended due to e-beam irradiation. If the e-beam 
is scanned across a vortex, the vortex's magnetic-field dis- 
tribution is distorted in the direction of the beam spot; 
i.e., the center of the distorted field distribution is dis- 
placed by 8r, which depends upon the distance of the 
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FIG. 1: (a) SQUID washer design; dotted line indicates grain 
boundary, (b) 5$ image of washer SQUID. The arrow indi- 
cates the direction d of one of the vortex signals. Tick labels 
are in units of /im. (c) Line scan 5&(x) along the dashed 
line shown in (b). The signal contrast C, i.e., the difference 
between maximum positive and negative signal from a single 
vortex, is indicated by the arrow. 



beam spot from the vortex. This displacement changes 
the amount of stray magnetic flux <I>(r) in the negative 
z direction that a vortex at position r couples into the 
SQUID hole. Hence, scanning across a vortex induces a 
negative (positive) flux change (5$ in the SQUID if the 
vortex is moved away from (towards) the SQUID hole. 
Figure IHb) shows an example of a 8$ image with vor- 
tices appearing as pairs of positive (bright) and negative 
(dark) signals. 

Fig- dc) shows a linescan 5$>(x) of the vortex signal 
along the dashed line shown in Fig. ^b). The difference 
between the maximum positive and negative signals from 
a single vortex defines the contrast C of a vortex in the 
(5$ image. The contrast C depends on the maximum vor- 
tex displacement Ar, which is independent of the vortex 
position, and on the gradient of the function $(r). Ne- 
glecting the (small) variation in |V$(r)| within Ar, we 
can write 



C = 2Ar |V$(r) 



(1) 



We define the direction d of a vortex signal 
[c.f. Fig. ^b)] as the direction of a unit vector in the 
(x, y)-plane pointing from the maximum negative to 
maximum positive signal of the vortex. Hence, d is par- 



allel to the gradient of <E>(r) and given as 
V$(r)/|V*(r)| . 



d 



(2) 



In the vortex-free case and for currents below the 
critical current of the Josephson junctions, the SQUID 
washer acts as a closed superconducting loop, which can 
carry a persistent supercurrent /, flowing clockwise or 
counterclockwise. The sheet-current density J(r) is di- 
vergence free; thus it can be written as 



J(r) = Iz x VG(r) , 



(3) 



where z is the unit vector in the direction perpendic- 
ular to the film plane and G(r) is a stream function. 
G(r) depends on the Pearl length, the film thickness and 
the washer geometry; it can be chosen such that G = 
on the outer edge and G — 1 on the inner edge of the 
SQUID[13- This situation corresponds to the "trapped- 
flux case" (zero applied field, magnetic flux trapped in 
the hole and slit, which is bridged at the edge) shown in 
the upper left plot of Fig. 1 in Ref. Htl 

While the function $(r) applies for the vortex state of 
the SQUID washer and the function G(r) applies for the 
vortex-free state, Clem and Brandt used an energy argu- 
ment to show that both functions are closely related |13| 
as 



G(r) 



$(r) 



(4) 



where $>Q=h/2e is the magnetic flux quantum. This 
equation holds as long as the current distribution in the 
vicinity of a vortex is not significantly altered by other 
vortices. This requirement is well fulfilled in our exper- 
iments since the applied fields are small (up to 40 fiT) 
and thus the distance between the vortices is a few /xm 
or more, well above A w 0.5 at T = 77 K. 
From equations one obtains 



Tf \ C 1 

Jir) = ^-2Ar 



(z x d) , 



(5) 



i.e., | J(r)\ oc C and J_Ld. The magnitude of the vortex- 
free sheet-current density J at the position r of an im- 
aged vortex is proportional to the contrast C of the LT- 
SEM signal given by that vortex, while the direction of 
J is perpendicular to the direction d of the vortex signal. 
These properties allow us to use the pinned vortices as 
local probes for the magnitude and direction of J. 

In order to check the validity of © we calculated the 
scalar stream function G(r) and the vortex-free sheet- 
current density J(r) for the geometry of our SQUID 
washers, usin g th e numerical simulation software pack- 
age 3D-MLSI [ljj. For the washer geometry in Fig.^a), 
G(r) and J(r) are shown in Fig. H{a) and (b), respec- 
tively. 

In Fig.[3]we compare experimental data obtained from 
four 5$ images at different cooling fields (upper row) 
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with the numerically calculated sheet-current distribu- 
tion J(r). In Fig. EJe-h) (middle row) every vortex 
imaged by LTSEM is represented by a white arrow at 
the position of the imaged vortex. The direction of this 
arrow is perpendicular to d, and the length is propor- 
tional to the signal contrast C. From Eq. (JSJ it follows 
that the arrows should indicate the direction and magni- 
tude of the normalized vortex-free sheet-current density 
J /I. For comparison, Fig. EJi~ni) (bottom row) shows 
the results from numerical simulations; i.e., the black ar- 
rows represent the direction and magnitude of J(r)/I. 
As shown in Fig. El the agreement between experimen- 
tal data (middle row) and numerical simulation (bottom 
row) is very good. The rms deviations in magnitude and 
direction of the sheet-current density are given in Tabled 
for four values of the cooling field. Since the proportion- 
ality constant between C and J /I contains Ar, which de- 
pends on experimental parameters such as e-beam power 
and sample temperature, we use Ar as a fitting param- 
eter for each <5$ image [c.f. Table H] to obtain the best 
agreement between experimental and numerical simula- 
tion data. There seems to be a trend towards smaller 
values of Ar with increasing B . The origin of this has 
not been clarified yet. 

The deviations between experimental and simulation 
data for J/I can easily be explained by the difficulty to 
extract the direction and magnitude of the vortex sig- 
nals exactly. Furthermore, vortices imaged close to the 
Josephson junctions yield a somewhat larger deviation 
in the direction of the sheet-current density as obtained 
from numerical simulations. During LTSEM imaging, 
shifts of the sample position up to a couple of /im in 
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FIG. 2: Numerical simulation for SQUID washer in the 
vortex-free case: (a) Scalar stream function G; G = 1 (white) 
at the slit edge; G = (black) at the outer edge, (b) Stream- 
lines, showing the current distribution in the SQUID washer. 



TABLE I: Rms deviation between experimentally determined 
and numerically calculated sheet-current density J(r) and 
values for maximum beam-induced displacement Ar. 



cooling field 
(MT) 


rms deviation 
of magnitude 


rms deviation 
of direction 


Ar 
(nm) 


5 


16% 


8.4° 


103 


10 


20% 


6.4° 


85 


15 


19% 


10.7° 


82 


40 


26% 


7.6° 


79 



both the x- and y-directions may occur. While such a 
shift along the x-direction can easily be corrected using 
the left and right (vertical) edges of the SQUID washer 
as a reference, a shift along the y-direction cannot be cor- 
rected by using the bottom washer edge as a reference, 
since we avoided directly hitting the Josephson junctions 
with the e-beam. In particular, for the largest cooling 
field B = 40 /xT, the analysis of the LTSEM data yields 
a significantly larger horizontal component for |J| as com- 
pared to the numerically calculated current distribution 
[c.f. Fig. Elf)]. Most likely, this deviation is due to a 
vertical shift which occurs during LTSEM imaging. Ex- 
cept for the signals close to the Josephson junctions we 
did not find a systematic deviation between experimental 
and numerical data. 

In conclusion, we have shown that the relation between 
the vortex-generated flux $ in a SQUID washer and the 
scalar stream function G allows us to use the vortices as 
local detectors for the sheet-current distribution J in the 
vortex-free case. This is a novel and more direct method 
of imaging circulating supercurrents that avoids compli- 
cated calculations. The experimental data obtained in 
magnetic fields up to 40 f/T are in excellent agreement 
with numerical calculations of J, confirming both the va- 
lidity of our model describing the generation of the LT- 
SEM vortex signals and the validity of the relationship 
derived by Clem and Brandt [l3j , even in the presence of 
many (up to 200) vortices in the SQUID washer. 
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Agreement (GIF) No. G-705-50. 14/01. D. D. gratefully 
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FIG. 3: Upper row: LTSEM 5$ images for different cooling fields Bo- Middle row: Magnitude and direction of sheet-current 
density calculated from the corresponding <5<E> images (white arrows) . Bottom row: Magnitude and direction of sheet-current 
density calculated by numerical simulation of J(r) (black arrows). Tick labels are in units of /im. 
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